Fluctuations and correlations of a driven tracer in a hard-core lattice gas 
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We consider a driven tracer particle (TP) in a bath of hard-core particles undergoing continuous 
exchanges with a reservoir. We develop an analytical framework which allows us to go beyond the 
standard force-velocity relation used for this minimal model of active microrheology and quantita- 
tively analyze, for any density of the bath particles, the fluctuations of the TP position and their 
correlations with the occupation number of the bath particles. We obtain an exact Einstein-type 
relation which links these fluctuations in the absence of a driving force and the bath particles density 
profiles in the linear driving regime. For the one-dimensional case we also provide an approximate 
but very accurate explicit expression for the variance of the TP position and show that it can be 
a non-monotoneous function of the bath particles density: counter-intuitively, an increase of the 
density may increase the dispersion of the TP position. We show that this non trivial behavior, 
which could in principle be observed in active microrheology experiments, is induced by subtle 
cross-correlations quantified by our approach. 

PACS numbers: 05.40.Jc, 05.40.Ca 



I. INTRODUCTION 

Active probe-based microrheology, which monitors the 
position of a tracer particle (TP) driven by an external 
force, has become a powerful experimental tool for the 
analysis of different systems in physics, chemistry and 
biology (see, e.g., Ref. [JJ). In the constant-forcing mi- 
crorheology setup, a microscopic particle embedded in 
a sample is actively manipulated by applying a known 
force. One measures then the response of the TP and po- 
tentially the microstructural deformations of the medium 
to learn about its microrheological properties. From the 
theoretical point of view, an unresolved issue, which is 
disregarded in the available continuous analytical ap- 
proaches, is to take into account explicitly the discrete- 
ness or "granularity" of the medium [JJ. This aspect 
becomes crucial when the probe and the medium parti- 
cles have comparable sizes, e.g. in colloidal suspensions. 
In this regime, both fluctuations of the probe that can 
not be described correctly if the medium is treated as a 
continuous bath and super-diffusion regimes have been 
reported [2H5J 

Within a broader context, modeling the response of 
a medium to a perturbation created by a tracer parti- 
cle, biased by an external force, is a ubiquitous problem 
in physics, which has been the subject of a large num- 
ber of theoretical works [7H5]. The resulting stochastic 
dynamics of the whole system is however a many-body 
problem which is difficult (or even impossible) to solve 
even in the simplest case when the particle-particle in- 
teractions are a mere hard-core. For that reason, in most 
approaches the microscopic structure of the bath is not 
taken into account explicitly, and the response functions 
are determined instead by using some effective bath dy- 
namics, modelled via Langevin or generalized Langevin 
[TO] equations (see also Refs. [ZHS])- While these ap- 
proaches are rather efficient, they cannot account for the 



detailed correlations between the tracer particle and the 
fluctuating density profiles of the bath particles. In this 
paper, we show that these cross correlations are of crucial 
importance and can lead to non trivial effects beyond the 
usual analysis of the force-velocity relation. We predict 
in particular a non monotonic behavior of the variance 
of the tracer position with the density of the bath, that 
in principle could be observed in active microrheology 
experiments. 

Our analysis relies on a model of driven tracer diffu- 
sion in a hard-core lattice gas, which appears as a min- 
imal model of active microrheology that explicitly takes 
into account the dynamics of a bath of discrete particles: 
A TP driven by an external force performs a biased lat- 
tice random walk in a bath of hard-core particles, which 
themselves perform symmetric random walks constrained 
by the condition of a single occupancy of each lattice site. 

At the theoretical level, the one dimensional version of 
this model is related to several well known models of out 
of equilibrium statistical physics. For example, in the ab- 
sence of external driving, it reduces to the single file diffu- 
sion problem [11) . In the absence of external driving and 
when additional adsorption/desorption processes with a 
reservoir of particles are considered, it corresponds to the 
so-called dynamical percolation [T21[T3]- Last, in the case 
of a constant external forcing experienced by all the par- 
ticles, it identifies with the asymmetric exclusion process, 
which has now become a paradigmatic model, both in the 
absence (see [H] for a recent review) or in the presence 
[T5] of adsorption/desorption processes. This model of 
driven tracer diffusion has been investigated both in the 
physical ,16"2fJ and in the mathematical (2TJ [22] litera- 
tures. Most of the obtained results, including proofs of 
the Einstein relation, are limited to the large time be- 
havior of the mean position (X tr ) of the tracer particle 
and the stationary density profiles of the bath (r) r ), where 
r\ T stands for the occupation number of the site r of the 



lattice, equal to 1 if occupied by a bath particle and 
otherwise. A general method to analyze the fluctuations 
of the tracer position (SXf r ) and the cross-correlation 
functions of the form (5X tr 8rj T ) is still lacking |23j . 

Here, for this minimal model of active microrheology 
we develop a theoretical framework allowing to analyti- 
cally determine these fluctuations in all regimes of bath 
particles density, and in principle for arbitrary dimen- 
sion. More precisely, our main results obtained within 
this formalism are (i) An exact Einstein-type relation sat- 
isfied by the cross-correlations in the absence of a driving 
force and the density profiles of the bath particles in the 
linear-response regime, (ii) An explicit approximate but 
accurate expression for the variance of the tracer position 
in a one-dimensional case valid for any density, (iii) A 
finding that the variance of the position of the tracer can 
be a non-monotonic function of the density of bath parti- 
cles, meaning that, counter-intuitively, an increase of the 
density of hard-core particles can increase the dispersion 
of the tracer position. We show that, in fact, subtle cor- 
relations between the tracer and the bath particles are 
responsible for this intriguing behavior. We anticipate 
that this striking effect could in principle be observed 
experimentally in the context of active microrheology. 



II. THE MODEL 

Consider a <i-dimensional hypercubic lattice of spac- 
ing a in contact with a reservoir of particles kept at a 
constant chemical potential (see figure [I]). Suppose next 
that the particles in the reservoir may adsorb onto va- 
cant lattice sites at a fixed rate f/r*. The adsorbed 
particles may move randomly along the lattice by hop- 
ping at a rate l/2dr* to any of 2d neighboring lattice 
sites, which process is constrained by a hard-core exclu- 
sion preventing multiple occupancy of any of the sites. 
The adsorbed particles may desorb from the lattice back 
to the reservoir at rate g/r*. The occupancy of lattice 
sites is described by the time-dependent Boolean variable 
r] r , which assumes two values, 1, if the site r is occupied 
by an adsorbed particle, and 0, otherwise. Note that the 
mean density of the bath particles, (rj T ), approaches as 
t — » oo a constant value p s = //(/ + g) but the number 
of particles on the lattice is not explicitly conserved in 
such a dynamics. 

At t = we introduce a TP, whose position at time t 
is denoted as R tr (t). The TP dynamics is different from 
that of the adsorbed particles in two aspects: first, it 
can not desorb from the lattice and second, it is subject 
to an external driving force E, which favors its jumps 
along the direction corresponding to the unit vector e! 
of the lattice. Physically, such a situation is realized in 
the context of active microrheology where the force on 
the TP is classically exerted by magnetic tweezers pQ. 

The TP dynamics is defined in the usual fashion: We 
suppose that the tracer, which occupies the site R tr at 
time t, waits an exponentially distributed time with mean 




FIG. 1: Model notations in the two-dimensional (2D) case. 



r, and then attempts to hop onto one of 2d neighboring 
sites, Utr + e„, where e v are 2d unit vectors of the hy- 
percubic lattice. The jump direction is chosen according 
to the probability p v , which obeys: 



p v = exp 



P 



(E-e„) 




exp 



P 



(E • e„) 



(1) 



where j3 is the inverse temperature and J^ denotes sum- 
mation over all possible orientations of the vector e M . 
The hop is instantaneously fulfilled if the target site is 
vacant at this moment of time; otherwise, i.e., if the tar- 
get site is occupied by any adsorbed particle, the jump 
is rejected and the tracer remains at its position. 



III. EVOLUTION EQUATIONS 

We focus on the projection of the TP position X tr = 
R tr • ei, which evolves in an infinitesimal time interval 
di according to 

X tr (t) + a, w.p. 7Ti =pi(l - r/i )dt 
X tr (t + dt) = { X tr (t) - a, w.p. 7r_i = p_i(l - T7_i)dt 

X tr (t),W.p. 1-7T1 -7T_1 

where w.p. stands for with probability. Averaging this 
equation, it is easy to find that the mean (X tr ) and the 
variance (SXf r ) of the TP position verify: 

~{X tr (t)) = ~{Pi(l - fe) -P-iO- - *-i)}, (2) 

at t 



and 



±(5X tr (t) 2 ) = y(l-p s ) 



(3a) 



T 

2a 



-{pi{ki - p s ) + p-i(k-i - p s )} (3b) 



{Pi9i -P-19-i} 



(3c) 



where 



k\ = (r/Rtr+x) and 9x = (^X tr 5rm. tr +x) (4) 



and 



fc±i = fc± ei and g±i = g± e 



(5) 



Note that the right-hand side of (3a) is a trivial mean- 
field term in which the effective jump frequency is given 
by (1 — Ps)/t. In turn, the two other contributions are 
nontrivial: the first one (3b) originates from the differ- 



ence between the actual density of bath particles in the 
TP neighborhood and their unperturbed value p s . The 
second one pc) involves the correlations between the TP 
position and the occupation numbers in front and behind 
the tracer. As we proceed to show, these two non trivial 
contributions, referred to in the following as the density 
and the cross-correlation contributions, are responsible 
for a non trivial behavior of the variance (5Xf r ). 

Now, in order to calculate (X tr ) and (6Xf r ), we have 
to determine the density profiles k±\ and the cross- 
correlation functions ]j±i at the sites adjacent to the TP, 
which requires, in turn, the evaluation of the density pro- 
files k\ and the cross-correlation functions g x for arbi- 
trary A. The evolution equations of these latter quan- 
tities can be deduced from the master equation (see ap- 
pendix [A} but these equations are not closed. Actually, 
they are coupled to the higher order correlations func- 
tions, so that one faces the problem of solving an infinite 
hierarchy of coupled equations. Here we resort to the 
simplest non-trivial closure of the hierarchy, by using the 
following decoupling approximation: 

(VWeJ « {VxiiVeJ (6) 

(SX tr r, xVe J » (r, x )(6X tr 5 Ve J + (r) eii )(6X tr 8rix)(7) 

for A 7^ {0, ±ei, ±e2 . . . , ie^}. The approximation in 
Eq. ffity has already been used in Rcfs. [18 20 and was 
shown to result in very good quantitative estimates of the 
TP mean position, while the one in ufh is new and, as we 
proceed to show, gives approximate expressions for the 
variance of the TP position which are in excellent agree- 
ment with numerical simulations. Note that these rela- 
tions are natural mean-field type approximations, since 
they can be regarded as expansions at the leading order 
in the fluctuation parameter Sy x = i] X — (r/x). 

Relying on these approximations, we finally obtain 
closed equations for k x and g x which hold for all A, ex- 
cept for A = {0,±ei,±e 2 . . . , ±e<f}: 



and 



2dr*d t k x = Lk x + 2df, 



2dr* 



(8) 



2dr*d t g x = Lg x + -' — er{pi(l - k ei )k x+ei 

T 



P-l(l - ^-eJ^A-eJ — ^2p^9 e ,y^k X 

2dT* 



2dT* 
T 



a{ Pl (l-k ei )-p^(l-k_ ei )}k x , (9) 



where L is the operator 

L = Y,A^^-2d(f + g), 



(10) 



V M is the discrete gradient defined by V M /(A) = /(A 

e ^) - fW, and 






(11) 



For A = e„ with v = {±1,±2, . . . , ±d}, i.e., the sites 
adjacent to the TP, we find 



2dr*d t k ev = {L + A v )k ev + 2df 



(12) 



and 



2dT*d t g e „ = {L + A v )g e 



2dr* 



a{pi(l - k ei )k Bv+ei 
. 2dr* „ , 

P-l(l - fc_ e jfce„-eij Pu9e u k 2ev 

T 

2dT* „ , 2dT* v-- „ _ , 

V-v9~^,k ev 2^ Pn9e^v^k ev 

<j{pi(l - fcei) -P-lO- - fc-ej}fce„- 

(13) 



T 
2dT* 



Note that Eqs.( 12 ) and ( 13 ) represent the boundary con- 
ditions for the general evolution equations dsb and pi, 
imposed on the sites in the immediate vicinity of the TP. 
Equations |8])-(13) together with Eqs.(2J) and pi) consti- 
tute a closed system of equations which suffices for com- 
putation of the variance (SXf r ) for any density and in 
arbitrary dimension. These equations also give access to 
the cross-correlation functions (5X tr 5r) X ). 



IV. EXACT EINSTEIN-TYPE RELATION 

Our starting point here is the usual Einstein relation 
p — /3D between the mobility p and the diffusion coef- 
ficient D, which has been shown to hold exactly for the 
system we study [55]. We now use the relations ([2| and 
(pi) in the long time limit to obtain the mobility and the 
diffusion coefficient in terms of the density profiles and 
the correlation functions: 

p = lim ^-{pi(l ~ fci) -p-i(l - fc-i)} (14) 



and 



D = 



lim 



2dT 



{Pl(l-fcl)+p_l(l-ft_l)} 



dT 



{P\9i ~P-i9-\\ 



(15) 



Using EqQ, as well as the trivial symmetry relations 
g-i{E = Q) = -g 1 (E = 0) and k ±l = p s ± aE + o{E) 
where a is a non vanishing constant, we finally get the 
exact Einstein-type relation of the form 

(3(6X tr S mtr+ei )(E = 0) = lim Q {m ^ } ~ Ps . (16) 



This equation relates a cross correlation function (be- 
tween the tracer and the bath) in the absence of field, to 
the linear response of the bath itself, and is compatible 
with the fluctuation dissipation relation described in [5]. 



V. EXPLICIT SOLUTION IN ID 

Lengthy but straightforward calculations lead to ana- 
lytical expressions of the density profiles k\ and cross- 
correlation functions g~\, allowing in particular to deter- 
mine the dispersion coefficient K tr = lim t _ J . 00 (<5X t 2 r )/2t 
(expressions given in appendix \Bj . Figure K2] shows that 
this expression of K tr is in excellent agreement with re- 
sults of numerical simulations for any density, which val- 
idates the decoupling approximation in Eqs. (|6| and tff\. 
We here discuss three important consequences stemming 
from this expression. 




words, the approximations (fsl) and (l7| are compatible 
with each other, which is a further validation of their 
applicability. 

Second, in the high density limit, p s — > 1, the disper- 
sion coefficient obeys 



K tr 



/ + v//(/ + 2) 



p^i2t/ +v //(/ + 2 ) + 2t*/- 



-(!-/>.)• (18) 



At leading order in 1 — p s , K tr is thus independent of 
the amplitude of the driving force. Quite unexpectedly, 
it turns out that the dependence on the driving force of 
the density and of the correlation contributions to the 
dispersion coefficient exactly compensate each other in 
this limit. 

Last, our analytical expression shows that, strikingly, 
K tr is a non-monotonic function of the density p s , for 
values of g smaller than some critical value g c , and has a 
marked maximum at a certain value of p s (see figure ^| . 
This result is a bit counter-intuitive, since one naturally 
expects that the dispersion will be maximal when there 
are no hard-core bath particles (i.e., when p a = 0). 
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FIG. 2: Dispersion coefficient as a function of the density p 3 . 
Analytical expression (solid lines) vs numerical simulations 
(symbols) for different biases (filled symbols : pi — 0.98 ; 
empty symbols : p\ = 0.8) and desorption rates g (•, o : 
0.15 ; A : 0.2 ;■ : 0.3 ; : 0.8). Inset : The different 
contributions to the dispersion coefficient, as defined in the 
main text, as a function of the density p s for g — 0.15 : trivial 
mean-field (dashed), density (dotted), and cross-correlation 
(dashdotted) contributions. The sum of these contributions 
is the dispersion coefficient (analytical expression represented 
by a solid line vs numerical simulations denoted by symbols). 

First, in the important limiting situation where no ex- 
ternal driving is applied to the system, we find that 



KtriE = 0) 



2t 



PsT' 



r U+g) i+^i+2(i+t*(i- Ps )/t)/(/+s) 



(17) 



Beyond its own interest, this expression together with 
the approximate expression for the velocity obtained in 
|18j in the linear driving regime using the decoupling ap- 
proximation (K3]) shows that the Einstein relation (exact 
for the system under study 22 ) is satisfied under the 
decoupling approximations, Eqs. ([6]) and ¥l\. In other 



FIG. 3: Cross-correlation functions g\ as a function of the 
distance to the tracer A for different values of the density p a 
and of the bias p\ . Analytical expressions (empty symbols) vs 
numerical simulations (filled symbols), (a) / = 0.033, g = 0.3 
; (b) and (c) / = 0.3, g = 0.3. 

To gain some insight into this non trivial behavior, we 
represent (see inset of fig [2]) three different contributions 
to the dispersion coefficient defined above (trivial mean- 
field, density, and cross-correlation contributions). As 
can be immediately seen, the non monotonic behavior of 
the variance originates from the non-monotonic behavior 
with p s of the cross-correlation function gi in front of the 
tracer. These cross-correlations functions g\, which pro- 
vide a quantitative measure of the correlations between 
the fluctuations of the tracer position and the occupa- 
tion numbers of the bath particles, are thus crucial in 
the problem. Actually, they display other non trivial 
and important behaviours (see fig. p| revealed by our 
analytical expressions. 

We first stress that the sign itself of these quantities 
is not trivial. While the analytical expressions allow us 
to show that, past the tracer, g~x is always negative and 
corresponds to an anti-correlation of the fluctuations, it 



is found that g\ is positive in the high density limit 
2p_icr(l — p 8 ) 



fji 



*-* 2 + [/ + y/f(2 + f)]r/i 
or in the absence of a driving force 



>0 



(19) 



9i(E = 0) 



all 



Ps 



-K tr (E = 0) > (20) 
a 



where K tr (E = 0) is given by Eq.(17), but can change 
sign otherwise. In other words, depending on the am- 
plitude of the driving force and the density of the bath 
particles, positive or negative correlations between the 
fluctuations of the tracer position and the occupation 
numbers of the bath particles can take place. In addi- 
tion, the monotonicity with the distance of these cross- 
correlation functions is also non-trivial - past (see fig. 
una)) and in front of (see fig. pTb)) the tracer - which 
further mirrors the complexity of these correlations. 



VI. CONCLUSION 

In conclusion, we have studied the dynamics of a driven 
tracer in a bath of hard-core particles in the presence of 
adsorption/desorption processes. For this minimal model 
of active microrheology, we have developed an analytical 
framework to go beyond the standard force- velocity rela- 
tion and quantitatively analyzed the fluctuations of the 
position of the tracer and their correlations with the oc- 
cupation number of the bath particles at any position, 
for any density of the bath particles. This approach al- 
lowed us to obtain an exact Einstein-type relation be- 
tween these correlations in the absence of a driving force 
and the density profiles of the bath particles in the linear 
driving regime. In the one-dimensional case we also pro- 
vided an approximate but accurate explicit expression of 
the variance of the position of the tracer and shown in 
particular that it can behave non-monotically with the 
density of bath particles: counter-intuitively, increasing 
the density of hard-core particles can significantly in- 
crease the dispersion of the tracer position. We have 
shown that subtle correlations between the tracer posi- 
tion and the occupation of the bath particles past and 
in front of the tracer are responsible for this non trivial 
behavior. Altogether, our results put forward cross cor- 
relations as observables to analyze active microrheology 
data and predict a striking non monotonic behavior of 
the variance of the position of the tracer with respect to 
the density of the bath, that could in principle be ob- 
served experimentally. In a broader context, Equations 
^ to (13) together with Eqs.Q and ^ provide a gen- 
eral framework to investigate important extensions that 
include the determination of explicit solutions in higher 
dimensions, in confined geometries and in non-stationary 
regimes. 
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Appendix A: Master equation 



Let rj = {?7r} denote the entire set of the occupation 
variables, which defines the instantaneous configuration 
of the adsorbed particles at the lattice at time moment 
t. Next, let P(R tr , 77; t) stand for the joint probability of 
finding at time t the TP at the site Rj r and all adsorbed 
particles in the configuration 77. Then, denoting as rf ,v 
a configuration obtained from 77 by the Kawasaki-type 
exchange of the occupation variables of two neighboring 
sites r and r + e„ , and as if a configuration obtained 
from the original 77 by the replacement 77 r — > 1 — rj T , 
which corresponds to the Glauber-type flip of the occu- 
pation variable due to the adsorption/desorption events, 
we have that the time evolution of the configuration prob- 
ability P(R tr ,?7;i) obeys the following master equation : 



2dT*d t P{R tr ,Tr,t) = 
d 

J2 Yl {p(Rtr,v r ^;t)-P(Rtr,mt)} 

M=l r/R tr - e M ,R tT . 

ft 
-(l-r ?Rtr+eM )P(R tr ,r7;t)} 



-2d 9 J2 {(l-Vr)P(Ktr,r, r ;t)-VrPC&tr,T,t)} 

r#R tr 

-2df Y, {llrP(Rtr,fi r ;t)~(l-Vr)P(Rtr,r,;t)}- 

(Al) 



r/R t , 



Note that the time evolution of the first moment 
(X tr {t)) (Eq. [2| can be obtained by multiplying both 
sides of Eq. (Al) by (R tr • ei) and summing over all 
possible configurations (R tr ,77). 



Eqs. [8] and 12 (respectively [9] a nd [l3| are obtained 
by multiplying both sides of Eq. ( |Al[ ) by 77(R tr + A) 
(respectively X tr r)(R t r + A)), summing over all possible 
configurations, and using the approximation R^ (respec- 
tively [7]). 



Appendix B: Explicit expression of K± 



2. Determination of g\ and g_i 



1. Determination of ki and fc_i 

According to Eq. [3] and defining A t; 

lim t _>oo<<«&>/2t, 



2 2 

Ktr = 2r ^ ~ ^ ~ 27 ^ 1 ^ 1 ~~ Ps -* 
+p_i(fc_i - p a )] (pigi -p-ig-i) 

T 



(Bl) 



Thus, the computation of k±\ and g±\ from Eqs. [8 12 
[9j and [l3| gives the expression of K tr . 

For one-dimensional lattices, a general solution of Eqs. 
[8] and [12] has the following form: 



k =k(X)-\ Ps+K+r,i forn>0 < f B 2) 

kn - k[X) -\ Ps +K_r2 forn<0. (B2) 



where 



'"? = 2A~ x \ A ^+ A -^ + 2 (f + 9) 



± 



^/(A + ^-i 



2(/ + .g) -4A!A_! 



(B3) 



while the amplitudes A± are given respectively by 



#4. = 



+ = P 



A_ = 



A x - A_! 

M - A_i 
'A-ifa-Ai 



(B4) 
(B5) 



Now, we are in position to obtain a system of two 
closed-form non-linear equations determining implicitly 
the unknown parameters A\ and A_\, which will allow 



us to compute k\ and fe_i. Substituting Eqs. B2 B4 and 
|B5| into Eq. [TT] we find 



Ai = 1 
A-! = 1 



2 P lT* 

T 
2 P -lT 



1 - Ps - Ps 



1 - p a - p. 



A x -A_ x 
A_i/n - A x 



A_i - Air 2 



(B6) 

(B7) 



Using Eqs. [9] and 13 and the expressions B2 of k n , the 
general equations satisfied by g become 

Ai(g n +i ~ 9n) + -A-iCfln-i - flw) - 2(/ + ff)flf„ 

+ ^a jpiA+r? ( 1 - p. - A+n - ^ J (n - 1) 

- P-xK+rl (l - Ps - K+n + ^J (rf 1 - 1) J = 

(B8) 
for n > 1, 

Ai(g n +i - g n ) + ^-i(? n -i - g n ) - 2(/ + g)5« 
+ ^a {pi A_r 2 " (l - p 8 - K +n - 9 -±\ (ra - 1) 

- p_i#_rj (l - p s - A+n + ^j (r 2 - J - 1) J = 

(B9) 
for n < — 1, 



/ 2r* 

A1.92 - .91 U-i + 2(/ + 5) + — pi(p s + A+r?) 

2r* _ 

+ P-i(ps + K + n)g-i 

T 

= -— opi(l - As - K +ri )( Ps + K+r{) 



2t* 



o"[pi(l-p« - i^+ri) 



p_i(l - p s - A_r 2 x )](p s + A+n), 



(BIO) 



and 



2r* 



A-iflf-2 - .9-1 A t + 2(/ + flf) + p_i(p s + A. 



2t* 

T 

2t* 

T 

2r* 



Pi(p s + K-r^Jgi 

op_i(l - p. - K-r^ 1 )^ + K_r^ 2 ) 

cr[pi(l -p«- A'+n) 



- p_ x (l - p s - A-_r 2 - x )](p s + A^ 1 ). 

The general solution is written 

1/ 



(Bll) 



(B12) 



i rW? for n > 0, 

/3r 2 l zrr nr 2 f° r n < 0- 

A x r 2 - A_tr 2 

where a and /? are constants to be determined, 



V = K+^n i Pl M - p s - A+n - ^ )( n - 1) 



- p-i [i- p a - K+n 



9-i 



(rr 1 -!) (B13) 



and 



It* 



W = K_—a \ Pl ( 1 - p s - K+n 



fli 



(r 2 - 1) 



P-i 



I- p s - K + n + -—^- ) (r 2 



(B14) 



Substituting Eq. B12 into Eqs. BIO and Bll on the 



one hand, and writing Eq. B12 for n = 1 and n = — 1 on 
the other hand, leads to a linear system of four equations 
satisfied by the four unknowns a, j3, ~g\ and g~i- It is 
finally found that: 



CE-BF , „ AF-CD 

■9i = ~TR 777^ and 9-1 = 



AE-BD 



AE-BD' 



(B15) 



where 



r Airi-A^ir/ 



A X T X 



- (a- X + 2(/ + g) + ^ Pl (p s + K + rl)\ (B16) 



D = 



£ 



and 



2r* , N A_ir~^ 

a Pl K_(r 2 - 1 ^-i_ 

r Air 2 - A_ir 2 

2t* 



+ Pifo + ^.rJ 1 ), 

r 



(B19) 



i4j + 2(/ + g) + ^P-iiPs + K_r^ 2 )\ (B20) 



It* 



A x v\ 



+ — p.iOo. + lf+ri), (B17) 



C 



It* 



ct(pi(1 -p s - if+ri) 



P-l(l - P. - ^-^2 1 ))(P« + K + r l) 



It* 

t 

2r* 



P-l 1 - Ps - 



crpi(l - p s - if+ri)(/9 s + A + r 2 ) 

pi(l - p., - A + ri)A+ri 

A ir f 



'2 



»~1 



Axn-A-irr 1 ' 

(B18) 



F 



2t* 



ff(pi(l 



ff+n) 



p_i(l - p. - ^-r 2 - 1 ))(p s + A^ 1 ) 

2t* 

op_i(l - p, - A-?^ 1 )^ + A-r^ 2 ) 

T 

2t* r 

cr (pi(l-p s - K + r 1 )K_r 2 



P-i 1 



A_\ A_ 
r"2 / r 2 



A_ x ro 



Air 2 - A_ir 2 



(B21) 



where A + , A_, Ai, and j4_i have been determined in 
the previous section. 
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